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The Hubbard model, which augments independent-electron band theory with a single parameter to 
describe electron-electron correlations, is widely regarded to be the ‘standard model’ of condensed 
matter physics. The model has been remarkably successful at addressing a range of correlation 
effects in solids, but beyond one dimension its solution is intractable. Much current research aims, 
therefore, at finding appropriate approximations to the Hubbard model phase diagram. Here we 
take the new approach of using ab initio electronic structure methods to design a material whose 
Hamiltonian is that of the single-band Hubbard model. Solution of the Hubbard model will then 
be available through measurement of the material’s properties. After identifying an appropriate 
crystal class and several appropriate chemistries, we use density functional theory and dynamical 
mean-field theory to screen for the desired electronic band structure and metal-insulator transition. 
We then explore the most promising candidates for structural stability and suitability for doping 
and propose specific materials for subsequent synthesis. Finally, we identify a regime - that should 
manifest in our bespoke material - in which the single-band Hubbard model on a triangular lattice 
exhibits exotic d- wave superconductivity. 


I. INTRODUCTION 


Independent-electron band theory - in which interac¬ 
tions between electrons are treated at a mean-field level - 
is a simple yet often remarkably effective method for cal¬ 
culating the behavior of electrons in solids. In its mod¬ 
ern density functional implementation, its combination 
of explicit description of crystal structure and chemistry 
with computational affordability has led to many suc¬ 
cesses ranging from early predictions of phase stability 
and lattice dynamics in technologically important ma¬ 
terials such as semiconducting silicon 1 , through descrip¬ 
tion of the structure and electronics of lattice defects 2 , 
to modern predictions of new physics and phenomena in 
complex materials such as nano- and heterostructures 3 . 

By construction, however, independent-electron the¬ 
ories do not include explicit correlations between indi¬ 
vidual electrons and so are often inappropriate for de¬ 
scribing or predicting the properties of so-called strongly- 
correlated materials, in which these interactions domi¬ 
nate the physics. Its classic failure is the inability to 
predict an insulating state for a half-filled band, which 
will always be metallic within the conventional band the¬ 
ory formalism. Here model Hamiltonians with explicit 
elect on-electron interaction terms have an advantage, at 
the expense of loss of computational tractability and 
chemical information. In particular, Hubbard introduced 
the modern form of the model that now bears his name in 
a series of papers beginning in 1963 4-6 specifically to “set 
up the simplest possible model containing the necessary 
ingredients”. (We note that the first hints of the Hubbard 
model can be found in the quantum chemistry literature 
in the 1950s 7 9 , with P. W. Anderson also suggesting an 
early variation 10 , and the standard form also indepen¬ 
dently proposed by Gutzwiller 11 and Kanamori 12 ). The 


Hubbard Hamiltonian is 

N 

H = —t ( C i,cr C j,CT + c j j(T c i,a) + U • (1) 

<z,j)cr i= 1 

Here the first term is the kinetic energy term of stan¬ 
dard band theory, with (i, j) indicating summation over 
all pairs of lattice sites (usually nearest neighbours), a 
indicating the spin degree of freedom, and t the hop¬ 
ping matrix element between states i and j. c\ and q 
are the usual creation and annihilation operators, and 
the number operator, rii = c]q. The second term intro¬ 
duces explicit electron-electron correlations through the 
Coulomb repulsion, of magnitude [/, between two elec¬ 
trons occupying the same lattice site. We illustrate both 
terms schematically in Fig. 1. 

For small values of U and large values of t (the weak- 
correlation limit), Eqn. 1 yields metallic solutions for par¬ 
tially filled bands and can successfully describe, for ex¬ 
ample, itinerant magnetism in transition metals. Impor¬ 
tantly, increasing the ratio oiU/t (the strong-correlation 
limit) causes a crossover to an insulating solution, al¬ 
lowing the description of, for example, the antiferromag¬ 
netic insulating state of transition metal oxides. Mod¬ 
ern studies can incorporate (non-self-consistent) informa¬ 
tion about crystal structure and chemistry via density 
functional calculations of the t and U parameters, and 
an increasing range of materials from Mott insulators 
to high-T c superconductors is currently addressed 13 16 . 
While it is not obvious why a single-band picture is rel¬ 
evant in the case of systems with degenerate d orbitals 
or strong hybridization with 5 or p states, several works 
show that the single-band picture does indeed success¬ 
fully describe low-energy charge and spin excitations in 
a realistic manner 17,18 . 

Despite the simplicity of Eqn. 1, the general Hubbard 
Model is in fact computationally intractable, and much 
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current research activity aims at obtaining the Hubbard 
model phase diagram. Exact solutions are possible in 
one dimension 19 , where a transition from a metal to a 
Mott insulator occurs at half filling when U > t 20 . Away 
from exactly half filling the solution is always metal¬ 
lic for either hole or electron doping. In two dimen¬ 
sions, numerically exact solutions can be found in the 
ground state only for either very small lattices of at 
most about twenty sites 21 , or on so-called ladder models 
consisting of coupled chains for around 100 sites using 
the density matrix renormalization group algorithm 22 . 
For arbitrary dimension and large lattices, however, ap¬ 
proximations must be used to investigate the phase di¬ 
agram. Analytical methods such as slave-boson mean 
field theory and the Gutzwiller approximation have been 
applied extensively 23-25 , and numerical techniques such 
as the Dynamical Cluster Approximation (DCA) and 
Dynamical Mean-Field Theory (DMFT) have had some 
success 14,26-28 . The rapid scaling of the computational 
resources required with the system size has hampered ef¬ 
forts in simulating large systems, however. Despite this 
profusion of computational approaches, consensus has yet 
to be reached on the exact nature of the phase diagram 
throughout the entire phase space, although diverse nu¬ 
merical methods start to agree on some properties of the 
2D Hubbard model 29,30 . Much less is known in three di¬ 
mensions, where substantial interest has been raised by 
the possibility of realizing the Hubbard model in optical 
lattices 31-33 and simulations have mainly been performed 
by quantum Monte Carlo simulations at relatively high 
temperatures of the order of the Neel temperature 34 39 . 

In this work we take a new approach to “solving” 
the single-band Hubbard model and use ab initio elec¬ 
tronic structure methods, at the band theory and dy¬ 
namical mean field theory levels, to design a real ma¬ 
terial whose Hamiltonian is exactly that of the single¬ 
band Hubbard model. Measurement of the properties of 
our bespoke material would then allow extraction of the 
Hubbard model phase diagram. Our approach is some¬ 
what reminiscent of experimental simulations of the Hub¬ 
bard model that have been performed using ultracold 
atoms 40 , which take the role of the the electrons in the 
conventional Hubbard model. The atoms are trapped 
in a potential created by interfering laser beams and 
since the lasers are highly controllable, a vast range of 
potentials can be created using different interference ef¬ 
fects. The simulated lattice model is therefore highly 
tunable both in terms of the lattice geometry and di¬ 
mensionality, and in the value of t/U. In particular, 
a metal-insulator transition, accompanied by a suppres¬ 
sion of doubly-occupied sites, a drop in compressibility 
and an excitation gap, has been demonstrated in a gas 
of cold fermionic atoms by tuning the strength of their 
mutual repulsion 41,42 , and short range magnetic corre¬ 
lations have been observed 43,44 . While proposals have 
been made to cool below the Neel temperature, however, 
(see, for example Ref. 37) to achieve long range magnetic 
oder and adiabatically prepare d-wave superconducting 


phases 45 , experimental progress has been slow and chal¬ 
lenging. A condensed matter realization of the Hubbard 
model, even if not as tunable as ultracold quantum gases, 
will thus be interesting as it will allow much lower tem¬ 
peratures to be reached. 


II. DESIGN OF A SINGLE-BAND HUBBARD 
MODEL MATERIAL 

A. Candidate crystal structures and chemistries 

Since we seek a material described by a single¬ 
band Hubbard model, we begin by engineering a non¬ 
degenerate d -manifold using crystal-field considerations. 
Crystal-field splittings with single non-degenerate d 
bands are obtained for a variety of ionic geometries, 
such as pentagonal bipyramidal, square antiprismatic, 
square planar, square pyramidal and trigonal bipyrami¬ 
dal (Fig. 2). In this work we choose the trigonal bipyra¬ 
midal coordination, in which the five-fold degeneracy is 
broken into two doubly-degenerate orbital levels (d xz , d yx 
and d xy ,d x 2 _ y 2 ), and a single orbital (d^), for its com¬ 
bination of large energetic separation of the single-band 
state and geometric simplicity. For trigonal bipyramidal 
coordination with shorter apical than in-plane ligand dis¬ 
tances, the isolated d 2 band is the highest energy level, 
and so it becomes half filled for a d 9 electron configura¬ 
tion. 

Such trigonal bipyramidal coordination is found in 
the well-studied hexagonal manganite structure adopted 
by YMnOs as well as the manganites of the smaller 
rare earth ions. The structure is formed from layers 
of corner-sharing MnOs trigonal bipyramids alternating 
with planes of Y ions. Figs. 3(a) and (b) show the high- 
symmetry PQ^/mmc variant that occurs at high tem¬ 
perature; below around 1000K, there is a phase transi¬ 
tion to a ferroelectric P 63 C 771 structure, which does not 
alter the coordination environment and hence does not 
change the crystal-field splitting. In Fig. 3(c) we show 
our calculated density of states for hexagonal YMnOs 
in the high symmetry PQ^/mmc reference structure, and 
with the local moments on the Mn ions ferromagnetically 
aligned. The 3d 4 Mn 34 ions have a high-spin configura¬ 
tion, with the up-spin d xzy d yx and d xy and d x 2_ y 2 or¬ 
bitals occupied, and the spin-up d z 2 forming the lowest 
energy state above the Fermi level. The exchange-split 
spin-down states begin at an energy of ^4 eV above the 
Fermi level. While the d 4 Mn 3+ configuration is magnetic 
and does not yield a half-filled single band, excluding it 
from consideration as our single-band Hubbard material, 
we see already that the d 2 orbitals form an isolated band 
with little hybridization with the oxygen ligands, con¬ 
firming our intuition that this crystal class is promising. 

Next we select ion combinations that have formal va¬ 
lence states that should yield half-filling of the isolated 
d 2 band, and will also likely allow for doping across the 
range from empty to filled band. As mentioned above, 
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to half-fill the d z 2 band, our B-site cation should have 
a d 9 configuration, suggesting Cu 2+ as the most promis¬ 
ing candidate. For oxides or sulphides, both of which 
have divalent anions, the A-site cations should be four 
valent suggesting Zr or Sn as possibilities. For fluorides, 
in which the anions are monovalent, the A-sites should 
also be monovalent, suggesting Li or Na. We take, there¬ 
fore the following as our initial list of trial compounds 
for further study: ZrCuOs, SnCuOs; ZrCuSs, SnCuSs; 
LiCuFs, NaCuFs. 


B. Computational details 

Our electronic structure calculations were performed 
within density-Functional theory using the Vienna ab 
initio Simulation Package (VASP) 46,47 . The electronic 
wavefunctions and density were expanded using a plane- 
wave basis set, and we used project or-augmented-wave 
(PAW) potentials 48 for core-valence separation. For the 
exchange-correlation potential we used the generalized 
gradient approximation (GGA) 49,50 and strong corre¬ 
lation effects were treated by means of the GGA+U 
scheme. Here we used the Dudarev method 51 in which 
a U e ff = U — J accounts for the on-site V (Coulomb 
replusion) and J (Hund’s rule exchange) on the metal d 
states. To choose an appropriate U e ff for hypothetical 
compounds (without spectroscopic data), we performed 
hybrid functional calculations to select an appropriate 
range by matching the DFT+U band gap with the hybrid 
functional result. In the hybrid functional calculations, 
an HSE functional consisting of 75% PBE and 25% ex¬ 
act Hartree-Fock exchange was used 52 . For our Cu-based 
compounds, this led us to choose a U e ff of 7 eV on the 
Cu -d states, which is similar to that used in previous 
studies on cuprate materials 53 . (The calculated YMn 03 
DOS shown in Fig. 3 (c) was also obtained using a U e ff 
of 7 eV on the Mn -d states, consistent with literature 
studies on YMn 03 .) A 10 x 10 x 4 Gamma-centred k- 
point mesh was used for Brillouin-zone integrations. The 
plane-wave cut-off was set to 550eV and in performing 
the structural optimizations, we allowed the ions to re¬ 
lax until the Hellmann-Feynman forces were less than 1 
meV/A -1 . 

To treat the electron-electron correlations in modeling 
the metal-insulator transition, we used the Dynamical 
Mean-Field Theory (DMFT) approach 54 . We first con¬ 
structed maximally localized Wanner functions for the 
transition-metal d bands of the Kohn-Sham Hamiltonian 
using the wannier90 code 55,56 . We then used these as 
the noninteracting part Hq of a Hubbard Hamiltonian 
H = Hq + Hi nt where Hi nt is the local electron-electron 
interaction comprising the intra-orbital Coulomb interac¬ 
tion, Udmft, and Hund’s rule coupling, J. Here we used 
J = 0 eV since we consider a single orbital, while we var¬ 
ied Udmft from 0 to 4 eV. We used a continuous-time 
hybridization expansion quantum Monte Carlo solver 57 
as implemented in the TRIQS 1.0 code 58 to calculate the 


local Green’s function within DMFT at a temperature 
T = 1 /(k B P) = 290K (/3 = 40 eV” 1 ). 


III. RESULTS 
A. Structural optimization 

First we perform an initial band structure screen¬ 
ing of our trial compounds within the high-temperature 
Pfts/mmc structure. Our calculated lattice constants for 
all six compounds are given in Table I. For the oxides 
and fluorides, the transition metal - apical ligand dis¬ 
tances are shorter than the corresponding in-plane dis¬ 
tances, which should yield the highest energy singly de¬ 
generate crystal field state required in our analysis. As 
expected from the larger atomic radius of Na, the vol¬ 
ume of NaCuF 3 is slightly greater than that of LiCuF 3 , 
although the apical Cu-F bonding distances are almost 
the same in both compounds. Between the Zr and Sn 
oxides, the volumes and Cu-0 distances are almost iden¬ 
tical, suggesting similar crystal-field splittings in the two 
cases. Interestingly, for the sulphides, the apical Cu-S 
bonding distances are longer than the in-plane distances. 
This should place the singly degenerate d 2 z band in the 
lowest energy position which will lead to its full applica¬ 
tion. Therefore we exclude the sulphides from our list of 
candidates for further consideration. 


B. Electronic properties at the density functional 
theory level 

Using our optimized structures, we next calculate the 
electronic properties of the candidate materials at the 
density functional theory level. Fig. 4 shows the calcu¬ 
lated electronic band structures and orbital-resolved den¬ 
sities of states for (a) ZrCu 03 , (b) SnCu 03 , (c) LiCuF 3 
and (d) NaCuF 3 . For both of the oxide compounds, the 
Fermi energy intersects a band that is composed primar¬ 
ily of Cn-d z 2 states, with some admixture from apical 
0-p z states. The Cu-d z 2 band is almost separated from 
the broad valence band composed of the other four ma¬ 
jority Cu -d bands and oxygen, but the separation is not 
complete due to its non-negligible bandwidth of 1.2 eV. 
For the tin case, we also find Sn -s states at Ep (Fig. 4 
(b)), which are undesirable for achieving a single-band 
model. 

We find a more promising situation in the fluorides 
(Figs. 4 (c) and (d)) both of which show a completely 
split-off half-filled Cu-d z 2 band. The reduced anion hy¬ 
bridization in these more ionic compounds gives a smaller 
bandwidth of ~ leV and makes the fluorides more single- 
band-like. (Note that our calculation is for a two- 
formula-unit unit cell, and the two Cu ions are slightly 
symmetry inequivalent so two Cu d 2 bands are shown.) 
The electronic structures of LiCuF 3 and NaCuF 3 are re¬ 
markably similar in the region surrounding the Fermi 
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level since there is no contribution from the A sites in 
either case in this region. The only clear difference is 
that the crystal-field gap between the split-off band and 
the other Cu -d states is slightly larger in NaCuF 3 
0.3eV) than in LiCuF 3 (~0.2 eV). 

In summary, we find from density functional calcula¬ 
tions that the cuprates of Zn and Sn, and the copper 
fluorides of Li and Na are promising candidate single¬ 
band Hubbard model materials when in the hexagonal 
manganite structure, with the fluorides having the more 
desirable properties. Next we continue with one of the 
fluoride representatives - LiCuF 3 - and investigate the 
effect of explicit electron correlations at the dynamical 
mean field theory level. 


C. Electronic structure at the dynamical mean 
field theory level 

To study the effect of explicit electron-electron corre¬ 
lations on the d 2 band, we next calculate the spectrum of 
LiCuF 3 using dynamical mean field theory (DMFT). We 
first construct the maximally-localized Wannier function 
(Fig. 5 (a)) on a Cu site for the energy window corre¬ 
sponding to the isolated Cu -d z 2 band using the Wan- 
nier90 code 56 . As expected the Wannier function has 
primarily d z 2 character with small p-symmetry “tails” 
extending onto the neighboring fluorine ions. 

Using this Wannier function as a basis, we next calcu¬ 
late the spectral function, A(uj) within DMFT, by ana¬ 
lytic continuation using the maximum entropy method 59 
Our results are shown in Fig. 5(b) for Udmft values 
ranging from 0 to 1.5 eV. We see that the system is 
metallic at V = 0 eV. As the interaction parameter V is 
increased, the spectral weight gradually shifts from the 
quasiparticle peak at uj = 0 eV to the upper and lower 
Hubbard bands. Eventually a classic Mott transition oc¬ 
curs when the electron interactions are strong enough 
to cause the quasiparticle peak to disappear completely 
and a gap to form. Consistent with our Tr G(/3/2) data 
(not shown), the system gaps above Udmft = 1.25 eV. 
We comment briefly on the apparent inconsistency in¬ 
troduced by our use of a non-zero U value in our origi¬ 
nal DFT calculations, combined with a second U in the 
DMFT study. To address this issue, we repeated our 
DFT calculations in the GGA (U = 0) limit and found 
the band structure and band widths to be very similar to 
the GGA+U case, with the DFT+U treatment acting as 
a scissors operator to isolate the d z 2 orbital. Therefore 
we expect only small quantitative changes in the Udmft 
value at the metal-insulator transition if the initial Wan¬ 
nier functions were constructed from results of the U = 0 
GGA calculation; construction of the Wannier functions 
would be complicated, however, by the poorly isolated 
bands. 


D. Structural stability 

As a first step to assessing the feasibility of synthesiz¬ 
ing our proposed compounds, we calculate the relative 
structural stabilities of a range of structural isomorphs 
of LiCuF 3 . We consider the following crystal structures 
all of which are known to exist for compounds with 
ABX 3 chemistry: high-temperature hexagonal mangan¬ 
ite (P6 3 /mmc), low-temperature ferroelectric hexagonal 
manganite (P6 3 cm), low-temperature nonpolar hexago¬ 
nal manganite (P3c), cubic perovskite (Pm3m), ilmenite 
(P3), the ‘NaCuF 3 ’ structure (which corresponds to the 
ambient experimental ground state, PI), the ‘CoGe0 3 ’ 
structure (C2/c), and two non-centrosymmetric struc¬ 
tural variants - rhombohedral (R3m) and tetragonal 
(P4/mmm) - found in BaTi0 3 . We apply both positive 
and negative hydrostatic pressure by scaling the lattice 
parameters, keeping the fractional internal coordinates 
of the ions fixed for each calculation to ensure that the 
original symmetry group is maintained. 

Fig. 6 shows the resulting calculated energies as a func¬ 
tion of volume for the various structural isomorphs of 
LiCuF 3 that were lowest in energy over the range of vol¬ 
umes studied. We find that the lowest energy structure 
is the ‘NaCuF 3 ’ structure, however for an 8 % increase in 
volume (corresponding to a 2 % increase in lattice con¬ 
stant) the hexagonal manganite structure is the ground 
state, suggesting that it might be possible to achieve it 
using tensile strain. At much larger lattice constants the 
C 2 /c CoGe0 3 structure becomes stable. 


E. Doping 

Finally, we explore routes to electron- and hole- dop¬ 
ing LiCuF 3 so that the entire range of band filling can 
be accessed. First we artificially add (remove) up to 
one electron per formula unit without changing the ion 
configuration, using a compensating background positive 
(negative) charge to prevent electrostatic divergence. As 
expected, the Fermi level shifts downwards (upwards) as 
electrons are removed (added), giving a completely un¬ 
filled Cu-d z 2 band when 1 e - per formula unit is removed, 
and a completely filled Cn-d z 2 band when 1 e - per for¬ 
mula unit is added. No qualitative change apart from 
the rigid shift in the Fermi level is found over the en¬ 
tire doping range; in particular the split-off single-band 
character is retained. 

An experimentally accessible possibility for hole dop¬ 
ing is to substitute fluorine (which forms F - , with one 
negative charge per atom) by oxygen (which forms O 2- , 
with two). Replacing one fluorine atom per unit cell 
with oxygen yields the chemical formula L^Ci^FsOi, 
which corresponds to removing half an electron per for¬ 
mula unit. We calculated the relative energies of the 
substitutional sites for the oxygen by replacing each of 
the two inequivalent fluorine ions with an oxygen ion, 
performing a full structural relaxation, then comparing 
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the final energies. We found that in-plane subsitution is 
lower in energy by ^0.3 eV per formula unit than apical 
substitution; this is favorable for our Hubbard material 
design as the in-plane ions have minimal hybridization 
with the split-off single band. The resulting density of 
states (DOS) for in-plane oxygen substitution is plotted 
in Fig. 7(b) (left panel), with the orbital projected O-p 
states shown in pink. Note that the electronic structure 
in the region of the Fermi level is remarkably similar to 
that obtained by removing electrons in our first doping 
calculations. 

Possible routes to electron doping LiCuF 3 are to re¬ 
place some Li with Be or to introduce F vacancies. Sub¬ 
stituting one Li ion with a Be ion introduces half an ex¬ 
tra electron per formula unit with the chemical formula 
LiBeCu 2 F 6 . The influence of Be substitution on the elec¬ 
tronic structure is shown in Fig. 7(b) (right panel), where 
we have plotted the Cu-d z 2 band separately from the 
other Cu -d states. While no Be states appear near the 
Fermi level, the substitution causes an orbital reordering 
whereby the region next to the Fermi level now com¬ 
prises contributions from the full d manifold. Removing 
one fluorine atom per unit cell to give the chemical for¬ 
mula Li 2 Cu 2 F 5 also corresponds to adding half an elec¬ 
tron per formula unit. We calculated the relative energy 
of vacancy sites by removing the fluorine ion from each 
of the two inequivalent fluorine in turn and performing 
a full internal relaxation. We found that apical fluorine 
vacancies are lower in energy than in-plane vacancies by 
~0.1 eV per formula unit; the resulting density of states 
for apical vacancies is plotted in Fig. 7(e) (right panel). 
Unlike the Be substitution case, the region surrounding 
the Fermi level retains its Cu-d z 2 character and behaves 
more like the rigid-band model of electron addition. 


F. Exotic superconductivity 

Since their discovery in the 1980’s 60 , the cuprate super¬ 
conductors have remained one of the most studied classes 
of materials, as well as one of the most elusive in reveal¬ 
ing a comprehensive understanding of their behavior 61,62 . 
In particular, while they have the highest critical tem¬ 
peratures of all known superconductors, the details of 
the pairing mechanism are still unknown. The relevant 
features for high-T c superconductivity in the cuprates ap¬ 
pear to be the quasi two-dimensional electronic structure, 
the parent antiferromagnetic compound with a range 
of accessible dopings, the spin-^ magnetic moment and 
the single band at the Fermi level. For example, the 
La 2 Cu 04 cuprate parent compound contains 3d 9 octa- 
hedrally coordinated Cu 2+ ions. A strong Jahn-Teller 
distortion lifts the energy of the two e g orbitals yield¬ 
ing a non-degenerate half-filled Cn-d x 2_ y 2 band which is, 
however, strongly hybridized with the p x and p y states 
on the neighboring oxygens 17 . The ground state is an¬ 
tiferromagnetic and insulating due to the combined ef¬ 
fects of exchange splitting and Mott electron repulsion. 


Recognition of these characteristics have led to attempts 
to design new materials that might enhance exotic su¬ 
perconductivity. A particularly clever suggestion was to 
use geometric engineering to crystal-field split the degen¬ 
erate e g orbitals rather than relying on the Jahn-Teller 
distortion 63 . A proposed candidate system was a layered 
superlattice of (LaNi0 3 /LaA10 3 ) in which the symmetry 
and strain of the superlattice would lift the degeneracy of 
the Ni e g states. While an intriguing idea, in experiment 
charge-transfer physics has been found to dominate the 
behavior and supress superconductivity. 

Our bespoke single-band Hubbard material, LiCuF 3 , 
exhibits all the necessary ingredients for exotic super¬ 
conductivity. In particular, the single band at the Fermi 
level is strongly two-dimensional, and Cu 2+ ions are 
spin-^ with antiferromagnetic interactions. The trian¬ 
gular lattice might lead to additional exotic behavior. 
Therefore, as a final analysis, we investigate whether 
a superconducting transition might occur in this com¬ 
pound and provide a prediction of its critical tempera¬ 
ture. For this we use a variant of the Dynamical Cluster 
Approximation (DCA) 64 , which is the cluster extension 
of the (DMFT) 54 . In this method, the original infinite 
lattice problem is transformed into solution of a finite¬ 
sized cluster with periodic boundary conditions, embed¬ 
ded in a self-consistent mean field. This transforma¬ 
tion is achieved via a coarse-graining procedure of the 
Green’s function, in which the Brillouin zone is divided 
into Nc patches on which the self-energy is assumed to 
be constant. The DC A is able to treat all short range 
correlations between the electrons in the cluster exactly, 
while long-range correlations outside the cluster are taken 
into account by the embedding self-consistent mean-field. 
We recently extended the method 65-67 , to allow deter¬ 
mination of the critical temperature for large clusters 
at a strong interaction strength and to treat continuous 
lattice self energies and vertex functions in momentum 
space, giving a more accurate description and allowing a 
more stringent investigation of the superconducting gap 
function in momentum space 67 . We use this new DCA+ 
method here. 

DC A has previously been used to investigate the 
pairing mechanism of the Cwoper pairs in the high-T c 
cuprates 68-74 . In addition, Chen et al. 75 recently inves¬ 
tigated the superconducting transition on a triangular 
lattice and discovered two divergent pairing susceptibil¬ 
ities, one with d xy symmetry and one with d x 2_ y 2 sym¬ 
metry, for a cluster of size 6 and U = 8.5 1. Furthermore, 
their phase diagram suggests that the critical tempera¬ 
ture rises as the doping is decreased, up to a possible 
T c « 0.06 1 (with t being the hopping parameter). This 
further motivates our study, as our LiCuF 3 compound 
has a half-filled band. 

We use the DCA+ method at V — St to solve a 16-site 
cluster, which has been used in earlier papers investi¬ 
gating the high-T c cuprate superconductors to allow for 
a direct comparison with these earlier calculations. In 
Fig. 8, we show our calculated leading eigenvalue A of 
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the matrix T y. If this leading eigenvalue A crosses 1, 
the pairing susceptibility y = Xo/(l — T y) diverges and 
a superconducting transition occurs. As in Ref. 75, we 
find two diverging superconducting modes, one with d xy 
and one with d x 2_ y 2 symmetry. In the inset, we show 
our finding that the leading eigenvalues are linear on a 
log-log plot, which indicates that the eigenvalue behaves 
as A oc a {T — Tc) 1 and is thus in the mean-field regime. 
This allows us to extrapolate to lower temperatures and 
find a T c « 0.08 by solving the equation A = 1. In or¬ 
der to understand why the d xy superconducting mode is 
favoured over the d x 2 _ y 2 , one can look at their corre¬ 
sponding gap functions <F, defined as the corresponding 
eigenvector, in momentum space. In Fig. 9 , we show the 
momentum-space structure of the gap function at the 
lowest matsubara frequency for respectively the d xy and 
the d x 2_ y 2 superconducting modes as well as a cut of the 
superconducting gaps along the Fermi surface. We see 


that the amplitude of the d xy mode is larger than that 
of the d x 2_ y 2 mode, which leads to a larger lambda and 
critical temperature. 

IV. SUMMARY 

We have identified a class of materials with the hexag¬ 
onal manganite structure whose low-energy behavior re¬ 
flects the single-band Hubbard Hamiltonian, and propose 
LiCuFs and NaCuFs as promising candidates. We used 
density functional theory, dynamical mean field theory 
and the dynamical cluster approximation to character¬ 
ize LiCuFs in detail and showed that it displays the 
expected Mott transition, as well as exotic d-wave su¬ 
perconductivity. We hope that our work will motivate 
experimental synthesis and characterization of this and 
related compounds 76,77 , as well as exploration of bespoke 
materials with other model Hamiltonians. 
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FIG. 1. Cartoon of a triangular lattice illustrating the two parameters that enter the Hubbard model: the inter-site hopping, 
£, and the onsite repulsion, U. 
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FIG. 2. Common coordination environments for transition-metal ions that yield at least one singly degenerate crystal-field level, 
(a) pentagonal bipyramidal, (b) square antiprismatic, (c) square planar, (d) square pyramidal and (e) trigonal bipyramidal. 



FIG. 3. (a) Side view of the hexagonal-manganite structure of YMnOs- Layers of corner-shared MnOs trigonal bipyramids 
are separated by layers of Y ions, (b) Top view of the hexagonal-manganite structure showing the triangular lattice formed by 
the Mn ions which are connected via shared oxygen ions, (c) Calculated density of states of hexagonal manganite YMn 03 in 
its high symmetry reference structure with the Mn ions ordered ferromagnetically. The orbitally-projected Mn-c^ states are 
shaded in purple and the Fermi energy, E^, is set to 0 eV. Majority (minority) spin states are shown on the positive (negative) 
y axis. 
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FIG. 4. Calculated band structures and total and orbital-projected densities of states (DOS) for (a) ZrCu 03 , (b) SnCu 03 , (c) 
LiCuF 3 and NaCuF 3 , all in the hexagonal manganite structure. The Fermi level is set to 0 eV and is marked by the dashed 
line in each of the plots. 




FIG. 5. (a) Maximally localized Wannier function in LiCuF 3 obtained from projection of the d z 2 band for one Cu site, (b) 
Calculated DMFT spectral function for a range of U values as a function of frequency uj. The Fermi energy is indicated to 0 
eV, and filled orbitals are shaded. 
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FIG. 6. Calculated total energy as a function of volume for the three calculated lowest-energy LiCuF 3 isomorphs. The NaCuF 3 
structure (blue shaded region and inset) is the most stable structure for small volumes, and the hexagonal manganite becomes 
the lowest energy structure for volumes greater than 63A 3 (green shaded region). Finally for very large volumes the CoGe0 3 
structure is stabilized (shown in orange inset). 
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FIG. 7. Left, a) Density of states of LiCuF 3 with half an electron removed per formula unit, b) Density of states for oxygen- 
substituted Li 2 Cu 2 F 50 i. Right, (c) Density of states of LiCuF 3 with half an electron added per formula unit, (d) Calculated 
DOS for Be-substituted LiCuF 3 . (e) Calculated DOS for LiCuF 3 with a fluorine vacancy. In all cases the Fermi level is set to 
OeV in each case, and orbital-projected DOSs are shown for the relevant atoms. 
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FIG. 8. The leading eigenvalue A of the matrix T x- A superconducting transition occurs when this leading eigenvalue A 
crosses 1 since the pairing susceptibility x — Xo/(l — T x) will diverge. Inset: A log-log plot of the leading eigenvalue versus 
temperature. The linearity of the plot suggests that A behaves as A oc a (T — T c ) 7 . This behaviour indicates that we are in the 
mean held regime. This allows us to extrapolate lambda and find a T c equal to oc 0.08. 



FIG. 9. The momentum space structure of the leading superconducting gap functions & dx y (left) and <f> d x 2 - y 2 (middle) at the 
lowest Matsubara frequency. By looking at their projections on the Fermi-surface (right), we can understand why the $ d ® 2 -i /2 
has a lower critical temperature. Its gap is simply not as big as the gap of & dx y in the negative region. Since the critical 
temperature is proportional to the size of the gap, one would expect a lower T c for <f> d x 2 - y 2 than for & dx y . 
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a 

(A) 

c 

(A) 

Cu-F (in plane) 

(A) 

Cu-F (apical) 

(A) 

ZrCuOs 

3.75 

9.33 

2.16 

1.81 

SnCu0 3 

3.78 

9.38 

2.18 

1.82 

LiCuF3 

3.47 

11.13 

2.00 

1.85 

NaCuF 3 

3.60 

11.60 

2.08 

1.84 

ZrCuS 3 

4.69 

10.26 

2.70 

2.22 

SnCuS 3 

4.72 

10.94 

2.71 

2.36 


TABLE I. Calculated structural parameters for hypothetical oxides, fluorides and sulphides in the centrosymmetric hexagonal 
manganite structure with the PQ^/mmc space group. The lattice parameters, a and c are given in the standard hexagonal 
setting with two formula units per unit cell. 




